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CN ■ Abstract 
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ly^ ' We proved the existence and uniqueness of a traveling wave so- 

0^ ■ lution to the thin film equation with a Navier slip condition at the 

liquid-solid interface. We obtain explicit lower and upper bounds for 
_ the solution and an absolute error estimate of approximation of a 

f"^ ■ solution to the thin films equation by the traveling-wave solution. 
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1 Introduction 

The degenerate parabolic equation 

ht + {h^h,,,), = (1.1) 

arises in description of the evolution of the height y = h{t, x) of a liquid film 
which spreads over a solid surface (y = 0) under the action of the surface 
tension and viscosity in lubrication approximation (see [U [H]). Lubrica- 
tion models have shown to be extremely useful approximations to the full 
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Navier-Stokes equations for investigation of the thin liquid films dynamics, 
including the motion and instabilities of their contact lines. For thicknesses 
in the range of a few micrometers and larger, the choice of the boundary 
condition at the solid substrate does not influence the eventual appearance 
of instabilities, such as formation of fingers at the three-phase contact line 
(see [5l[in|). For other applications, such as for the dewetting of nano-scale 
thin polymer film on a hydrophobic substrate the boundary condition at the 
substrate appears to have crucial impact on the dynamics and morphology 
of the film. 

The exponent n G M+ is related to the condition imposed at the liquid- 
solid interface, for example, n = 3 for no-slip condition, and n G (0, 3) for 
slip condition in the form 

y^ = ^/i"-2(t;^)^ aty = 0. (1.2) 

Here, v^ is the horizontal component of the velocity field, /x is a non-negative 
slip parameter, and /i /i"~^ is the weighted slip length. Distinguished are the 
cases /i = and n = 2. The first one corresponds to the assumption of a no- 
slip condition, the second one to the assumption of a Navier slip condition at 
the liquid-solid interface. The wetted region {h > 0} is unknown, hence the 
system is simulated as a free-boundary problem, where the free boundary 
being given by d{h > 0}, i.e. the triple junctions where liquid, solid and 
air meet. 

The main difficulty in studying equation (II. ip is its singular behaviour 
for /i = 0. The mathematical study of equation Hl.U was initiated by 
F. Bernis, A. Friedman [2j. They showed the positivity property of solu- 
tions to (|l.ll) and proved the existence of nonnegative generalized solutions 
of initial-boundary problem with an arbitrary nonnegative initial function 
from H^. More regular (strong or entropy) solutions have been constructed 
in (HE]. One outstanding question is whether zeros develop in finite time, 
starting with a regular initial data. What is known is that with periodic 
boundary conditions, for n ^ 3.5 this does not occur [HE], while for n < 3/2 
the solution develops zeroes in a finite time [6]. One way of looking at the 
problem (II. 3p has been to study similarity solutions to (11. ip in the form 
h{x,t) = t~"H{xt~^), where na + 4/3 = 1 (see |4j). In the paper [3], the 
authors proved also existence dipole solutions and found their asymptotic 
behaviour. We note that the solutions of such type do not exist in the case 
n ^ 2, however, there exists a traveling wave solution (see [8j). 

In the present paper, we concentrate on a traveling wave solution to 



sjj) s,(t) + w(t) 



Figure 1: Sketch of the cross-section of a dewetting film after rupture, 
showing the expanding dewetted residual layer (also called a "hole" or "dry 
spot") in the middle and the adjacent "dewetting ridges" moving into the 
surrounding undisturbed uniform film |12| . 



(II -ip at n = 2, namely, we consider the following problem with a regular 
initial function: 




(1.3) 

;a; = at X = Si. 

System (11.31) describes the growth of dewetted regions in the film. Fluid 
transported out of the growing dry regions collects in a ridge profile which 
advances into the undisturbed fluid (see Figure 1). Under ideal conditions, 
it could be imagined that dry spots could grow indefinitely large. By con- 
servation of mass, the growing holes would shift fiuid into the ever-growing 
rims. In our situation, large length scale to limit the sizes of these struc- 
tures is absent, and we might expect the motion and growth of the ridges 
to approach scale-invariant self-similar form. At the same time, the ridge 
profiles have a pronounced asymmetry (see [T2]). 

In the problem at hand, x = Si{t) is the position of the former moving 
interface, i. e. the contact line, while the position of the latter interface will 
give an effective measure of the width of the ridge, x = S2(t) = Si(t) +w{t). 
The ridge is assumed to be moving forward, si(t) > 0, corresponding to 
an expanding hole. The arbitrary positive parameter 9 corresponds to the 
contact angle of liquid-solid interface. Thus, we can control a dewetted 
region in the film by the contact line and obtain asymmetry profiles of so- 
lutions. As the paper [I2j has shown that the axisymmetric profile can 
be analyzed within a one-dimensional thin-film model. The authors found 



matched asymptotic expansion, speed and structure of the profile, in par- 
ticular, they obtained that 

h{x, t) ^ A si^/2(s2 - x)^/2 at X = sa, (1.4) 

where the asymptotic constant A = 2(2/3)^/^. 

Hereinafter, we assume that the contact line moves with a constant 
velocity (v) and the width of the ridge (w) is a constant. As in [I2], we are 
going to look for a solution to l\1.3h in the following form 

h{x, t) = h{^), where C, = x — vt, si = vt, S2 = vt + w . 

We can remove v and 9 from the resulted problem by rescaling appropriately. 

As a result, we obtain the following problem for the traveling- wave 

^(r7)v^'"(r7) = 1, ^iv) > 0, 

V?(0) = ^{d) = 0, (^'(0) = 1, 0^r]<:d. ^ ' 

Boatto et al. [8] reduced this problem to the problem of finding a co- 
dimension one orbit of a second-order ODE system connecting equilibria. 
Hence generically solutions will exist but only for isolated values of the free 
parameter d. The parameter d was found in [I2] by integration Hl.bL and 

d=l/2. (1.6) 

Our paper is organized as follows. In Section 2 we prove the existence 
and uniqueness of a traveling wave solution to the problem (11.51) (Theo- 
rem [T]) . Lower and upper bounds for the traveling wave solution are con- 
tained in Section 3 (Theorem [2]). We note that the bounds assert that 
the constant A of p.4p must be from the interval [4v^/3, 4-\/6/3] (Corol- 
lary 3.1). In Section 4 we find an absolute error estimate of approximation 
of a solution to (II. 3p by the traveling- wave solution (Theorem [3l). 

2 Existence of the traveling wave solution 

Below we prove the existence and uniqueness traveling wave solution to the 
problem (II. Sp . Our proof is based on some modification of the proof of the 
existence and uniqueness dipole solutions from [3]. 



Theorem 1. There exists a unique solution ip{rj) to the problem ( [J.5IJ such 
that (fi{ri) e C^{0, d) n C^[Q, d] and ^{rj) > for < rj < d. 

First, we prove the following auxiliary lemma: 

Lemma 2.1. Assume that ip G C^{0,d)nC^[0,d], ip{0) = ip{d) = 0, v?'(0) = 
1, V? > and (p"'{ri) ^ m {0,d). Then (p has a unique maximum and 
ip'{d) < 0. 

Proof. Since (p"'{ri) ^ we have that (p'{ri) is convex, (p"{ri) is increasing. 
By Rolle's theorem (p'{r]) has at least one zero in (0,(i), and (p'{r]) has no 
more than two zeroes in (0, d) by convexity. Let di,d2 G (0, c?) : (p'{di) = 
'p'{d2) = 0. Then, by Rolle's theorem, there exists d^ G ((^1,(^2) : '^''ids) = 
whence ip{d) ^0. In view of yj > in (0, d) and (p{d) = 0, we obtain d2 = d 
and (p'{d) = 0. This proves that f'^rj) has exactly one zero in (0, c?) and 
hence ip{r]) has a unique maximum. Now 'p'id) < follows easily. D 

Proof of Theorem\^ Green's function. We define a Green's function G(r7, t) 

by 

G"'{r], t) = 6{r]-t), ^ r] f^ d, O^t^^d, 

G(0, t) = G{d, 0) = G"(0, t) = 0, O^t^d, 
where d = 1/2. By explicit computation, we find that 



f2.11 



0<G(nt)-l 2(t-rf)Vif0^r]^t^rf, 

(,^U[r],t)-<^ 2(^t-d)W~d{r]-tfii0^t^r]^d, ^^■^> 



whence 



G{ri,t)dt=^ri'^{d-ri){l-ri), G{7],t) dt = \7f{d - 7]f (2.3) 



7) 

if < r/ < (i, and 

^(r],^) ^Ct2(t/-t), |G"(?7,t)KCt(rf-t) V77, tG [0,(i]. (2.4) 

Approximating problems. For each positive integer /c we consider the prob- 
lem 

Vkiv) = ^k^ for < r] < rf, 

MO) = Md) = i /fc(0) = 1, (2.5) 

Mv) e c3[o, rf], Mv) >0ioi0^ri^d. 

5 



Consider the closed convex set 

S = {ve C[0, d]:v^ 1/k in (0, d)} 
and the nonlinear operator $/j defined by 

d 

^kviv) = I + 2r/(rf -v)+ [ Giv, t)v-\t) dt 



where G{rj, t) is from (12.21) . The operator $fc mapping S into S is continuous. 
Moreover, ^k{S) is (for each k) a bounded subset of C^[0,(i] and hence a 
relatively compact subset of 5*. By Schauder's fixed-point theorem, there 
exists (fk G ^k{S) such that $fcV^fc = V^fc- This is the desired solution of the 
problem (12. 5p . Note that (pk satisfies 



d 

MV) = l + Md - V) + JG{r],t)^';!{t) dt, 



d 

cp',{7l) = 2{d-2^) + lG'{7i,tM'{t)dt. 





(2.6) 



In view of Lemma EH] (applied to (pk — ^/k), there exists a unique point rrik 
in which the maximum of ipk is attained. Therefore, 

iPkiv) / in (0,mfc) and ipkij]) \ in {mk,d). 
Estimates. Since G ^ we get 

Mv) - -k = Hd -V) + I G{v, t)ipl\t) dt ^ ipl\r^) J G{r], t) dt 



= Iv^d -v){l- v)Vk\v) > lv\d - vfVk'iv), 
whence 

^kiv) ^ ylv^^'^id-ri) ^ CT]{d-T]f^'^ if < r/ < m^, 
and 



d 

I 





Mv) - I = 2^(^ - ^) + / G{^, t)if,\t) dt ^ ^,\r]) J G{r], t) dt 

V 

= lv\d-7]Yip^\r]), 



whence 

(fikiv) > ylvid - Tlf^'^ if rrik <r] < d. 

Hence, 

Mv)^Cr]{d-r]f^ \/rie{0,d). (2.7) 

Next we deduce from the differential equation that, for all t] G (0, d), 

^{d-^)^';!{^) = ^{d-^)^l\r^) ^ C{d-r^)-^l\ (2.8) 

where the right-hand side is an integrable function. 

Passing to the limit. From (12.61) . (12.41) and (12.81) it follows that (pk is bounded 
in C^[0, c?]. Therefore, there exists a subsequence, again denoted by v^fc, 
and a function ip such that ipk^'V uniformly on [0, d] as A; ^ oo. Thus 
ip{Q) = ip{d) = and, by (12. 7p . </? > in {0,d). Hence, for each compact 
subset / of (0, d) we have 

and, by (12. 8p and Lebesgue's dominated convergence theorem, 

r]{d-r]M'{r])^ r]{d - r])^-\r]) in L\0,d). (2.9) 

Since ip'l! -^ ip'" in the distribution sense, it follows that 

fc— >oo 

^"'{7^)=^-\r])m{0,d), (2.10) 

i.e. ip satisfies the differential equation. Moreover, from (12. 9p . (12.10p . (12. 6p 
and (12. 4p we deduce that ipk~^f in C^[0,c?] as /c ^ oo, and hence ip also 
satisfies ip'{0) = 1. This completes the proof of the existence. 

Uniqueness. Let ipi and ip2 be two solutions of the problem ( ll.Sp and set 
V = (pi — ip2. Since v'" = (pi^^ — Lp2^ and the function v i— ^ v~^ is decreasing 
we deduce that vv'" ^ 0. Since vv"' ^ and f(0) = v{d) = v'{0) = 0, we 
conclude that v = 0. This completes the proof of Theorem [H D 



3 Lower and upper bounds for the traveling 
wave solution 

Integrating (11.51) with respect to r/, we arrive at the following problem 

^{r])^"{r]) = i(/(r7))2 + 7]-d, ip{v) > 0, ^g_^^ 



ip{0) = 0, /(O) = 1, ^r] ^d, 

where d is from p.6l) . Analyzing the behaviour of a solution to f l3.ip . we 
find explicit lower and upper bounds for the solution. 

Theorem 2. Let (p{ri) be a solution from Theoreml^ Then the following 
estimates are valid 

for all T] G [0, rf] (see Figure 2). 

Corollary 3.1. In particular, from ^3.^) it follows that 

¥^i'7^A-'2 - xf/^ ^ hix, t) ^ l^if |z^(s2 - xf' (3.3) 

for all X G [si, S2]. 

Lemma 3.1. The function i^oiv) = Aor]{d — r])^^'^ (A^ > 0) satisfies the 
inequalities 

if^" ^ liif')^ + 7]-dy7]e[o,d] if Ald^ ^ 5/3, (3.4) 

(^(^" ^ \{'^'f + r]-d^T]e [0,ci] if Ald^ ^8/3. (3.5) 

Proof. Indeed, the function ^Qirj) satisfies the equation 

where 

fir,) := lAlid -r,)(v- ^^^) {v ' ^^^^) ^ V., G [0, rf]. 

From 

f{ri) ^ri-d^bri^-Mri-Ad^ + jj^OyriE [0, d], 

D = Ad' + 20d' - f = f {Aid' - I) ^ if Aid"" ^ 5/3 

we obtain (13.41) . In a similar way, we obtain (13. 5p for /(r^) ^ rj — dWrj G [0, d] 
if Aid' ^ 8/3. D 
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Figure 2: Lower and upper bounds for the traveling wave solution 

Lemma 3.2. The function ipramiv) = Air]{d - 7]^^ {Ai > 0) if Afd'^ ^ 8/9 
is a lower bound for the solution ipirf) of ^3.1\) : 

Vminii]) ^ V^(^) («• e- Vminiv) " fijl) < 0) Vt] G [O, d]. 

Proof. "Contraction Principle". Let us define v{rj) := ipminiv) ~ fiv)- Sup- 
pose that there exists a point rjo G [0, d] such that ^(r/o) > then 770 is a 
point of maximum for v{ri), i.e. v'^tiq) = -v^ V^mm(^o) = v'iVo) = M and 
v"{t]o) < 0. From ([SID and ((331) we deduce that 



^"(v) ^ k f %!f^ 



(Vvf 



iv-d) 



whence 



fM>^=:M^ i^-M' +"»-''). ■ 



<0 



<0 



Using M = Ai{d- T]oy/\d - 5r/o/2), we find 

lM^ + 7]o-d={d~r]o)[lAl{d-5r]o/2f-l] ^OVt/qG [0,d] if A^rf^ ^ 8/9. 

Thus we obtain a contradiction with our assumption, which proves the 
assertion of Lemma I3.2[ D 

Lemma 3.3. The function ipn,e.^{r]) = A2r]{d-r]f/'^ {A2 > 0) if A'^d'^ ^ 8/3 
is an upper bound for the solution ipirf) of ^3.1\) : 

^iv) ^ V5max(^) («• e. ipir]) - V2max(^) ^ 0) Vr] G [0,d]. 

Proof. "Contraction Principle". Let us define f (r/) := ^9(77) — ipma^iv)- Sup- 
pose that there exists a point 770 G [0, d] such that f (r/o) > then r/o is 
a point of maximum for v{ri), i.e. f'(?7o) = -v^ V'max(^o) = V^'(^o) and 
v"{rio) < 0. Moreover, from v^max(^) — 'iAi{d — r])~^^'^[5r]/4: — d] we find that 
<^max(^) ^ for all 7/ G [0,0.4] and (fi'^^^iri) > for all r] G (0.4, rf]. From 
(13. ip and ( 13. 5p we deduce that 

r" il"' ^^7l ("J "^2 , „ _ - =^ ¥'(r]0)¥^"(r]0) - V'max(r/0Xax(^0) > 0, 
V^maxV^max ^ 2 V^^ max J "r '/ " 

(3.6) 
whence 

^^(^o)<^"(^o) - y^max(^o) <^max(^o) , = ViVoW {m) + V'L^iVo) v{Vol - 

This contradicts to our assumption if t^q G [0,0.4]. 

Now, let r/o G (0.4, d]. In this case, if (p"{rio) ^ 0, and we obtain a 
contradiction immediately from (13.61) . If ^"{tjo) > then we rewrite (13. 6p 
in equivalent form: 

^f^^^" Y W/ VA - d ^ (^'^^°) - ^n.ax(^0))" > 0, (3.7) 

whence 

i^<^^(^o)-ymax(^o))''^ = iv"im) +y'Lxivo))y 'im) 

^0 >0 <0 

and we arrive at a contradiction. Thus, Lemma 13731 is proved. D 

As a result of Lemmata 13.21 and 13.31 we obtain lower (more exact in 
comparison with (12.71) ) and upper bounds for the solution of (13.11) . and 
consequently for the solution of (II. Sp . 
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4 Absolute error estimate of approximation of 
a solution by a traveling-wave solution 

The next theorem contains an absolute error estimate of approximation of 
a solution (e.g., a generalized solution) by a traveling- wave solution. 

Theorem 3. Let h{x, t) be a solution and h{^) be the traveling-wave solution 
to the problem lll.3\) . Then the following estimates hold 

sup \h{x,t) - MO I < fe{s2 - si)i/2 tf \h,{s2)\ > e, 

X&[SI,S2] 

sup \h{x, t) - Mo - 2esl%2 - Sl)l/2| ^ ^Q^g^ _ ^^)l/2 -f |;^^(3^)| ^ Q^ 

xe[si,s2] 

(4.1) 
where ^ = x — vt, Si = vt and S2 = Si + w. 

Proof of Theorem H We make the following change of variables in Ill.Sfi 

h{x, t) 1-^ /(^, t), where ^ = x — vt. 

As a result, we obtain the following problem 

/*-v/5 + (/Vcce)« = 0, 

/(O, t) = f{w, t) = 0, /2/€?? = at e = 0, (4.2) 

f^{0,t)=e>0. 

Multiplying (I4.2li ) by —f(({C,,t) and integrating with respect to ^, we get 

W WW w 

H I /K^' i) ^e = -V / /,/,, d^ + j iPf,,,)d^, rfe = -i / |(/|) d^+ 



w 



+ 1 ilVmhh d^ ™^™^ lie^ - f^{w, t)) - / /V|,, d^, 




whence 



w w 

ki I f!i^^t)d^ + J ffl^d^ = lie' - fl{w,t)) ^ 



w 

lfl{U)d^^v{e'-fl{w,t)). (4.3) 



A 

dt 
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Integrating (14.30 with respect to time, we find 

t 
m^MWn,) ^ ll/«(e,0)||i.(o,.)+vy"(^^-/|(^,t))rft (4.4) 



From ( 1441 ) it follows that 

\\f(,m^ </ IIMe,o)|li.(o,.) if|/,(«;,t)|>^, 

From this, by virtue of uniqueness of the traveling-wave solution h{^) (see 
Theorem HD, /(^,0) = h{^) we deduce from (l45ll that 

ii/5(e,t) - /i^(oiiL2(o,^«) ^ 2\\h^{o\W{^,^) if i/e(«^,^)i > ^, 

ll/c(e,t) - /ic(OIU^(o,^«) ^ 2||/z5(OIU^(o,^«) + 20v^ if |/^(Ti;,t)| ^ e. 

o 

Therefore taking into account the embedding H^{^,w) C C[0, ty], we find 
that 

sup |/(e,t) - m\ ^ 2w^'^h^mL^io,^) if \h{w.t)\ > 0, 
?e[o,-io] 

sup |/(e,t) - m\ ^ 2w'/\\\h^mLHo,^.) + OWt) if \f^iw,t)\ ^ e. 

(4.6) 
Since /ig(0 = ^v'riiv) i'^iv) is from Theorem [T]) and (p{ri) has a unique 
maximum in (0, d), due to (13. 2p we conclude that 

||/^5(OIU^(o,«,) ^ §0. (4.7) 

Thus, we from (14.61) and (14.7P arrive at 

sup |/(e,t) - MOI ^ fow'/' if |/5(^,t)| > e, 

sup |/(e,t) - m - 29w'/'^/Vt\ ^ few'/' if mw,t)\ ^ 9, 

which completes the proof of Theorem [3l D 

Acknowledgement. Author would like to thank to Andreas Miinch for his valuable 
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05-1000008-7921. 
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